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ABSTRACT 

A number of strong infrared forbidden lines have been observed in several 
evolved Wolf-Rayet star winds, and these are important for deriving metal abun- 
dances and testing stellar evolution models. In addition, because these optically thin 
lines form at large radius in the wind, their resolved profiles carry an imprint of the 
asymptotic structure of the wind flow. This work presents model forbidden line pro- 
file shapes formed in axisymmetric winds. It is well-known that an optically thin 
emission line formed in a spherical wind expanding at constant velocity yields a flat- 
topped emission profile shape. Simulated forbidden lines are produced for a model 
stellar wind with an axisymmetric density distribution that treats the latitudinal ion- 
ization self-consistently and examines the influence of the ion stage on the profile 
shape. The resulting line profiles are symmetric about line centre. Within a given 
atomic species, profile shapes can vary between centrally peaked, doubly peaked, 
and approximately flat-topped in appearance depending on the ion stage (relative to 
the dominant ion) and viewing inclination. Although application to Wolf-Rayet star 
winds is emphasized, the concepts are also relevant to other classes of hot stars such 
as luminous blue variables and Be/B[e] stars. 

Key words: stars: emission line, Be - stars: mass-loss - stars: winds, outflows - 
stars: Wolf-Rayet 



1 INTRODUCTION 

Outstanding questions about the three-dimensional structure 
of stellar winds continue to plague our understanding of mass- 
loss processes from stars. Although much has been accom- 
plished toward understanding the gross properties of stellar 
winds in terms of theories for spherically symmetric mass- 
loss (a description of current understanding along with co- 
pious references in this regard can be found in Lamers & 
Cassinelli 1999), a complete physical description of two and 
three dimensional effects remains elusive. There have been no- 
table advances in a number of areas, with a brief and incom- 
plete listing to include such studies as axisymmetric models of 
magnetized thermally driven winds (Sakurai 1985; Washimi 
& Shibata 1993; Sean & Balick 2004), axisymmetric line- 
driven winds with rotation (e.g., Bjorkman & Cassinelli 1993; 
Owocki, Cranmer, & Gayley 1996; Petrenz & Puis 2000), 
axisymmetric magnetized line-driven winds (Poe, Friend, & 
Cassinelli 1989; Brown et al. 2004; Ud-Doula & Owocki 
2002; Townsend & Owocki 2005), axisymmetric dust-driven 



winds (Asida & Tuchman 1995; Dorfi & Hofner 1996; 
Soker 2000), and time-dependent spherical or time-averaged 
spherical models (Feldmeier 1995; Dessart & Owocki 2002; 
Runacres & Owocki 2002, 2005). The challenge has been 
that to include all of the relevant physics could require multi- 
dimensional radiative (magneto-)hydrodynamics. Despite the 
difficulties involved, semi-analytic and detailed numerical 
models have been inspired by a large volume of observations 
that exhibit behavior that the steady-state spherical models 
are unable to explain. The examples are numerous, such as 
the triple rings of the SN1987A remnant (Crotts, Kunkel, & 
McCarthy 1989; Burrows et al. 1995; Sugerman et al. 2005), 
the fact that the majority of planetary nebulae and luminous 
blue variables are aspherical (Balick 1987; Nota et al. 1995; 
Davies, Oudmaijer, & Vink 2005), the existence of discs in the 
Be and B[e] stars (see the review of Be stars by Porter & Riv- 
inius 2003; see the recent proceedings on B[e] stars by Kraus 
& Miroshnichenko 2006), discrete absorption components and 
related phenomena seen in hot star wind lines (Howarth & 
Prinja 1989; Massa et al. 1995; Brown et al. 1995; Lepine, 
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Moffat, & Henriksen 1996), and the inference of apparently 
ubiquitous dumpiness in stellar winds (Hillier 1991; Bouret, 
Lanz, & Hillier 2005; Fullerton, Massa, & Prinja 2006). 

The goal of this contribution is to consider diagnostics 
of aspherical mass-loss as revealed by forbidden emission line 
profiles. The formation of most forbidden lines occurs over a 
large-scale volume of the wind. Consequently, the emission 
profile shapes that form in the hypersonic winds of early-type 
stars are imprinted with the properties of the asymptotic wind 
structure. This is most certainly not a new concept. Dating 
back to Beals (1929), it was appreciated that flat-topped emis- 
sion profile shapes result for optically thin lines that form in 
a constant velocity and radially expanding flow. In a colliding 
wind binary, Stevens & Howarth (1999), Luhrs (1997), Hill 
et al. (2000), and Hill, Moffat, & St-Louis (2002) have shown 
how time-dependent deviations from flat-top shaped profiles 
can be used as a tracer of the geometry of the colliding winds. 

In this paper a parametric axisymmetric density distribu- 
tion is adopted as a framework for modeling forbidden emis- 
sion profiles as a function of the wind asymmetry and viewing 
inclination. In addition, we point out that the density function 
alone is not enough to completely determine emission pro- 
file shapes. The emissivity depends on the ionic species and 
its latitude dependence that is is coupled to the overall den- 
sity distribution. So our model profiles also include a treat- 
ment of the influence of latitudinal ionization effects. Interest- 
ingly, for a given stellar wind model, different lines can show 
different, even opposite, profile trends, as for example some 
lines show centrally-peaked emission, whereas others display 
double-peaked profiles. A number of simplifying assumptions 
are employed to focus on the general trends. The wind ex- 
pansion is assumed constant for large radius and also in lati- 
tude. Clumping is ignored. How clumping may vary with lati- 
tude and radius in an axisymmetric hot star wind has not been 
well-studied, and so we choose to avoid making ad hoc as- 
sumptions. However, the basic equations presented here can 
be modified to include velocity and clumping effects to study 
their impact on forbidden line profile shapes. 

Section 2 presents a brief overview of forbidden line for- 
mation in spherically symmetric stellar winds. Then Section 
3 contains a description of the axisymmetric model, and re- 
sults for model forbidden line profiles are given. Following 
is a summary and discussion of applications in Section 4. 
There are two appendices to conclude the paper. Appendix lAl 
presents an analytic solution for line profile shapes in the case 
of a pole-on viewing perspective of an axisymmetric wind, and 
AppendixlBlspecificallv shows that the two-level atom approx- 
imation that will be adopted throughout the paper applies to 
the particular and interesting case of Nelll for Wolf-Rayet star 
conditions. 



2 FORMATION OF FORBIDDEN LINES IN 
SPHERICAL WINDS 

Our presentation of the theory of forbidden line emission gen- 
erally follows the approach and to some extent the notation of 
Osterbrock (1989) regarding the atomic physics and Barlow, 
Roche, & Aitken (1988) in terms of stellar wind application. 



2.1 The Atomic Physics 

Consider a fine structure transition for an ion species assuming 
a two-level atom only, denoting ' 1 ' for the lower level and '2' 
for the upper level. The level population density (number per 
unit volume) m can be modified by collisional and radiative 
excitations, and the level population density ri2 can be mod- 
ified by collisional de-excitations and spontaneous radiative 
decays. Highly non-LTE conditions exist in hot stellar winds, 
and so the two-level atom representation can be an excellent 
approximation in many cases. However, it can be the case for 
some atoms that there are more than two levels that can be 
important for the population that determines the emission for 
the line of interest. In the case of the Wolf-Rayet (WR) stars, 
the lines of [Nell] 12.8 ( 2 P 1/2 - 2 -P3/2) and [Nelll] 15.56 
ixm ( 3 Pi - 3 P 2 ) are of considerable interest for determining 
the abundance of neon. Although, the ground state of Nell is a 
doublet and so should be well-described by the two-level atom 
approximation, Nelll has five low-lying levels. The particular 
line at 15.56 /im happens to involve the two lowest levels of 
the five, and we will show later that the two-level approxima- 
tion is pretty good for conditions relevant to WR stars. Thus 
the theory that follows assumes a two-level atom throughout. 

It is commonly the case that radiative excitations from 1 
to 2 are ignored compared to the collisional term. The radiative 
excitation rate (number per volume per second) is m -B12 J v , 
with B\2 the upward Einstein transition probability, and J v 
the mean intensity at the frequency of the line transition. The 
collisional rate is m n c 512, where n c is the electron number 
density and qi2 is the collisional excitation rate in volume per 
second. A ratio of these rates shows that the radiative excita- 
tion rate is indeed much smaller than the collisional rate. For 
the wind case, J v will decrease with radius in the wind as the 
dilution factor that scales as r~ 2 , and this cancels with the 
electron density in the denominator that will also scale in the 
same way. The dominant factor that remains is a ratio n c /n c .o, 
where n c is the 'critical density' and n c ,o (to be defined in 
the following section) is a scale factor for the electron number 
density. The critical density is defined here as 

A21 n\ 
n c = . (1) 

921 

In effect, the critical density represents a transition between 
high density regimes with (n c 2> n c ) for the line emissivity is 
linear in density versus low density regimes with (n e <C n c ) 
for which the emissivity scales with the square of density. 

The volume emissivity j (erg s _1 cm -3 sr _1 ) is given by 



3 = — h U21 A21 n 2 . 

4.7T 



(2) 



To determine 712, equilibrium conditions for the two levels are 
imposed as governed by the following rate expression: 



m n e <ji2 = n 2 n c g 2 i + n 2 A21, 
where 

<?21 91 3 



with 



<12 



P = hu 2 i/kT e , 



(3) 



(4) 



(5) 
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for V21 the frequency of the line transition, T c the electron 
temperature, and g the statistical weight of the level, with g = 
2 J + 1. The collisional volume rate is given by 



8.629 x 10" 6 fii 



<?21 



1/2 



92 



(6) 



where fii2 is the collision strength. 

We define n^E to be the number density of element E in 
ion stage i, and we approximate its value as 



ni,E « m + ri2. 



(7) 



We make the following definitions: the ion fraction for this 
element in this ion stage is 



«i,E . 



the abundance of the element relative to all nucleons is 

A - nE 
-4e — — , 

n N 



(8) 



(9) 



for tin the number density of nucleons; and the ratio of the 
nucleons to electrons is denoted 



7c = 



nn 



(10) 



With these relations, the emissivity in equation can be re- 
expressed in terms of n e instead of n 2 . First, equation is 
used to eliminate ni in equation 0, the result of which can be 
solved for 712 alone, giving 



"2 



m,B n c gi2 



n e <?12 + n c q 2 i + A 21 



(11) 



Factoring out n c qi 2 , using the definition of the critical density 
from equation Q in conjunction with the preceding relations, 
we arrive at 



n 2 = Qi,E Ae 7o n c 



1 + il e e + il^l e P 



32 n c 



(12) 



At high densities for which n c S> n c , one has that j oc A21 n e , 
whereas at low densities with n c <C n c , the scaling is j oc 
A21 nt/n c oc gi2 r? c . 

In Barlow et al. (1988), the quantity 71 is defined as the 
ionic fraction, being the number density ratio of a particular 
ion of a particular atom to all ions in the gas. In our notation, 
one has that 71 = Qi,E Ae- It is desirable to derive the ionic 
fraction to compare against stellar evolution predictions. The 
ionic fraction can be inferred from the observed total flux of 
forbidden line emission, provided that the source distance is 
known. If lines from different atomic species are available, line 
ratios can be used to derive abundances (or possibly limits to 
abundances) even if the source distance is not well-known. 



2.2 Spherical Winds 

It is instructive to briefly review the formation of forbidden 
lines in a spherical wind. For this purpose it is useful first to 
develop a sense of scale, to determine where the line emission 
is formed in the wind flow. 

Consider a spherical wind with mass-loss rate M and ra- 
dial speed v(r) for radius r. The mass density will be 




log r/R, 

Figure 1. A logarithmic plot of the radial gradient of the line lumi- 
nosity, dLi/dr oc r 2 j(r), versus normalized radius in the wind. 
The curve shows three main components. The first is the 'Accel- 
eration' component at far left where v(r) changes from some low 
value at the wind base to Vaa. Middle is the 'Flat' portion, where 
r 2 j oc constant. Right is the 'Linear' part where r 2 j oc r~~ 2 . The 
vertical scale is shifted so that the 'Flat' portion is at zero. The differ- 
ent breaks from the flat segment to the linear portions represent line 
transitions with different critical densities (see text). 



Psph = 



M 



4-7T r 2 v(r) 



(13) 



The number density will be n sp h = tin + n c = p sp h/m, 
with m = firriH the average particle mass in the wind, and 
the mean molecular weight. The electron number density 
is n c = Paph/ /J-efriH, for /i c the mean molecular weight 
per free electron. Similarly, the nucleon number density is 
tin = Psph/ Unitih, for /iN the mean molecular weight per 
free nucleon. In a hot star wind that is dominated by ionized 
hydrogen, the relations become m = mn/2, fi c = = 1, 
n c = and so "f c = 1. This would be appropriate for an 
O star, for example. Strong forbidden lines are seen in the IR 
spectra of WR stars, and a WR wind is often dominated by 
twice ionized helium at small radii but possibly once ionized 
helium at large radii (e.g., Fig. 3 of Dessart et al. 2000 shows 
Hell is dominant for their WC8 model of WR 135). Assum- 
ing that helium is once ionized, the molecular weights become 
/i = 2, /i c = 4, and = 4, and 7 C = 1. 

The region of line formation is roughly set by the ra- 
dius at which the electron density and the critical density are 
equal. Among hot stars the mass loss rates are around 10 -10 
Mq yr^ 1 and larger, with terminal speeds generally around 
1000 km s _1 or greater from spectral type B and earlier (e.g., 
Lamers & Cassinelli 1999). We define a number density scale 
factor to be 



M 



n 



4-7T R 2 Voa M m H ' 



(14) 
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with values from 10 7 cm -3 for B star winds soaring to 10 14 
cm -3 in the case of dense WR winds. The critical density for a 
forbidden line depends on the electron number density, so we 
define a related scaling factor for the electron density: 



(15) 



Typical stellar wind forbidden lines observed in WR winds 
have critical densities ranging from 10 4 to 10 7 cm -3 (e.g., see 
tab. 5 of Ignace et al. 2001). Ignoring the weaker winds with 
quite small M values, the radius in the wind where n c = n c 
occurs is expected to be very large. We can find this radius r c 
from n c = (n e ,o) (-R*/'"c) = n c giving 



n c ,o 

■ 



(16) 



with typical values of 10 2 — 10 4 stellar radii. 

A goal of studying forbidden lines in stellar winds is to 
derive information about metal abundances. The observable 
is the total line flux that is related to the line luminosity as 
Fi = Li /A-kD 2 , where D is the distance to the source. Being 
an optically thin line, the total line luminosity depends on the 
atomic rates associated with the line transition and on the wind 
parameters, including the ionic fraction. With the atomic rates 
known, and the wind properties such as mass-loss rate and ter- 
minal speed provided, the observed line emission can be used 
to find the ionic fraction. 

The total luminosity of the optically thin line emission is 
given by a volume integration involving the emissivity, viz 



Li = 4tt IT F, = 4tt 



j(r)dr. 



(17) 



Defining ui = ((71/(72) e , then substituting equation 1121 into 
equation J2j, and this again into equation 1171 . one obtains 



Li = Lo 



where 



dr/R* 



1 + u + u (r/r c ) 2 ' 



Lo = 4-7T Qi,E Ae 7e h V21 A21 n e ,o R* 



(18) 



(19) 



Using a change of variable with u = R,/r, the integral 
can be re-expressed as 



Li — Ln 



da 



l Q (1 + 0j)u 2 + UJU 2 ' 

where it c = R,/r c . This has an analytic solution: 



Li = io 



+ Uj) LOU 2 



tan 



1 /l +w 



(20) 



(21) 



With uj a factor of order unity, and the expectation that u c <C 
1, the total luminosity generated for this forbidden line is 



Lo 



tt/2 



V( 1 + uj ) uj 



nc 4 1/2 « 3/2 
— oc /i 21 n c , 



(22) 



which is in effect the same as equation (13) of Smith & Houck 
(2005), except that those authors allow for a constant clump- 
ing factor (see our discussion of clumping in sect. 13. 3> . From 



equation J22> one can easily backsubstitute to obtain the prod- 
uct of the ion fraction and metal abundance in terms of the 
measured line luminosity (for known distance): 



Qi,E Ae 



(1+w 



1 1 /2,„l/2 rl l/2 



Li 



2iy' 2 7 C h u 2 i A21 n 3/ 2 7? 3 



(23) 



The reader may object that our derivation is lacking an 
important component of the wind, namely the inner radii 
where the wind accelerates. We have so far tacitly assumed 
a constant expansion wind at v = Voo to derive the above 
result. At the inner wind, the flow accelerates from some near- 
hydrostatic configuration at subsonic speeds to the hypersonic 
speeds of order 10 3 km s . For a typical forbidden line and 
a fairly dense O or WR star wind, the value of the critical ra- 
dius is so large, that the contribution to the line flux from the 
very dense inner and accelerating wind region is negligible. 
To illustrate this, Fig. Q shows the contribution function for 
the line luminosity. Plotted is dLi jdr oc jr 2 representing the 
line emission from a spherical shell as a function of the ra- 
dius. A wind velocity law with the form v(r) = (1 — b/r) 
is assumed (with b defined so that v = vo some initial wind 
speed at radius r = i?*), such that no = 10 13 cm~ 3 , and 
n{R*) x 10 15 cm -3 . 

At small radii, the curve is sharply dropping because of 
the rapidly changing velocity law. This portion is labelled the 
Acceleration' branch. Interior to the critical radius (where 
n c ^> n c ), the contribution function is constant, and this 
branch is labelled as 'Flat'. At large radius beyond the criti- 
cal radius, the emissivity drops as the square of density, where 
r 2 j oc n e , and so is labelled 'Linear'. The different breaks 
from the flat portion are representative of forbidden lines with 
different critical densities, here ranging from 10 4 cm -3 (break 
at largest radius) to 10 8 cm -3 (break at smallest radius). 

Although the Acceleration branch shows quite a strong 
contribution function, it is relegated to such a narrow range of 
radii as compared to the location of the critical radius, that it 
can safely be ignored for computing the emission line flux and 
profile shape. What remains is a constant velocity spherical ex- 
pansion flow that forms a flat-topped emission profile, which is 
the standard theoretical result for this case of an optically thin 
line (e.g., Mihalas 1978; Lamers & Cassinelli 1999). Note that 
the Flat and Linear portions contribute approximately equally 
to the total line luminosity. 



3 AXISYMMETRIC WINDS 
3.1 Density Structure 

With a good picture of the forbidden line formation in spheri- 
cal winds in mind, the case of axisymmetric flows is now ad- 
dressed. To model such winds, the following parametric forms 
for the wind density are adopted. We assume the wind is top- 
bottom symmetric in what follows. Then if the equator has 
higher density than the pole, we use 



n(r,e*) =np(r) [1 + (G - 1) sin m 0,] . 



(24) 



where the star is taken to have spherical coordinates 
(r, #*,(/>*), n p (r) is the purely radius-dependent number den- 
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Figure 2. A plot of the effective opening angle A8 of a polar or 
equatorial wind plotted as a function of the power exponent m. 



sity distribution along the polar axis, and G 



q /n p > 1 is 



the equator-to-pole density contrast. The latitude dependence 
of the density is modeled as a sine function of the co-latitude 
6* to some power m. This density function is motivated by ex- 
amples of equatorial density enhancements inferred for some 
classes or particular stars, such as B[e] stars, or WR 134 that 
will be discussed later. In the case of polar enhanced wind 
flows, G = n p /n cq ^ 1, and we use 



n(r, 8 t ) = n cq (r) [1 + (G - 1) cos" 



(25) 



models that are inspired by systems like r\ Car and other Lu- 
minous Blue Variables that show bipolar morphologies. These 
density functions are not specifically related to any particular 
model, but are chosen as a convenient way of characterizing 
axisymmetric winds. 

To characterize the degree of asymmetry in the wind, an 
opening angle AO refers to the latitude interval between ei- 
ther the pole or equator and the latitude at which n(r, 9*) = 
0.5 (n p + n eq ). This works out in value to be the same for both 
polar or equatorially enhanced winds, as given by 

A8 = cos" 1 (2~ 1/m ) . (26) 

Values of A8 are plotted in Fig.[2]against the logarithm of the 
power exponent m. Within this parametrization, well-confined 
polar or equatorial flows require quite large values of m. For 
example, in order to achieve AO ~ 5°, a value of m ~ 10 2 is 
required. 

With forbidden line emission formed over a large-scale 
volume, we shall ignore the inner wind acceleration based on 
the arguments of the previous section, and treat all densities to 
vary as r~ 2 . To compare the spherical outflow case with the 
axisymmetric models, the following normalization is used: 



n sph (r) 



n(r, 8,) dfi t , 



(27) 



Figure 3. A plot of the latitude-dependent density normalized to the 
density that a spherical wind would have against stellar co-latitude. 
The curves are for polar enhanced (solid) and equatorially enhanced 
(dashed) winds for values of G = 3 and m = 3, or G = 10 and 
m = 33. 



where /i* = cos 8* . For a given value of G and m, this equa- 
tion determines the ratio H cq = n Bp h{r) /n p (r) for a dense 
equatorial wind or H p — n sp h(r)/n cq (Y) for a dense polar 
wind. 

To illustrate, example density profiles are plotted in Fig.|3| 
for different values of G and m. The density is normalized 
to the spherical equivalent, namely a spherical wind with the 
same mass-flux as the axisymmetric one. The density is plot- 
ted as a function of stellar co-latitude in cos and so the 
equator is at left, and the pole at right. The solid curves are for 
the polar enhanced winds, and the dashed for denser equatorial 
winds. A brief comment on the dense polar winds versus the 
dense equatorial winds. When m and G are constant, n p will 
be larger for a dense polar wind than the corresponding value 
of n eq for a dense equatorial wind. Although the opening an- 
gle A8 will be the same, the solid angle associated with the 
respective dense components is not equal, being smaller in the 
case of a polar flow. 



3.2 Ionization Structure 

For the ionization balance in the wind, steady-state condi- 
tions are assumed. Only radiative recombination and photo- 
ionization are considered. The ion stage i ranges from to 7, 
where '0' denotes a neutral atom, and 'J' complete ionization. 
Since interest is focussed on forbidden lines of trace metals, 
the electron number density is assumed not to be impacted by 
the ionization balance of the metals. For example in O stars, 
hydrogen can be safely considered as completely ionized with 
n c = hh for all radii and latitudes. For a WR wind, singly 
or doubly ionized helium would dominate the electron num- 
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ber density (although ionization of C and O is relevant for the 
carbon-rich WR stars). Mainly modest distortions of the wind 
from spherical are considered, and not dense equatorial discs 
or polar jets. The simplification is that with n c given, the ion- 
ization balance for different metal species is independent. 

Under these conditions, ionization equilibrium will be 
governed by equations of the form 



= +a i+ i(r c ) n c m+i - Qi(T c ) n c m 



+ I -^-a v (i — X)ni-\dv 
I hv 



a v {i) n\ dv, 

hv 



(28) 



where the ion stages are all for the same element so that the 
subscript 'E' has been dropped, ai is a recombination coeffi- 
cient for ion i, J v is a mean intensity, v\ is the photo-ionization 
edge frequency for ion i, and ov(i) is the photo-ionization 
cross-section for ion i. In general, the equation directly cou- 
ples three adjacent ion stages; implicitly all of the stages are 
coupled through the fact that he = J^. ni. To complete the 
set of equations, one has that for the equation involving i — 0, 
there is no recombination to a lower ion stage, and for the case 
i = I, there is no photo-ionization to a higher ion stage. 

These expressions can be simplified considerably, by as- 
suming the diffuse radiation is negligible. The application of 
the theory in this paper emphasizes WR winds, and it is well- 
known these winds are so dense that the continuum forms in 
the wind itself, and a classic static photosphere is not observed. 
As such, the stellar radius is ill-defined, being wavelength de- 
pendent. Still we adopt a core-halo approach, in which the con- 
tinuous opacity is thin beyond the radius where the continuum 
forms, because this is a good approximation for ascertaining 
the ionization balance at large radii where the forbidden line 
emission is generated. 

To determine the ionization balance from a set of expres- 
sions based on equation 1281 . it is useful to define the following 
dimensionless parameter: 



Si = 



n c ,o( 



where 



) ai+i 



H I //- 1 [1 + (G - 1) cos 11 
"" " [1 + (G — 1) sin" 



fie 



(29) 



(30) 



for the appropriate dense polar or dense equatorial wind case. 
For convenience, we introduce <5_ 1 = 1 . The factor T\ pertains 
to the photo-ionization rate, as given by 



hv 



(31) 



where the factor of 1 /4 comes from the dilution factor in the 
large radius limit. It is also seen that the latitude dependence 
for a wind with an axisymmetric density distribution enters 
into the ionization balance via the recombination rate terms. 
The ion fraction Qi,E will be given by 



Qi,E(0.) = 



(32) 



This form also holds for the Saha equation, such as found in 
Mihalas (1978, eq. 5.17). However, we stress again that equa- 
tion 1321 for the solution of the ionization balance specifically 
refers to trace metals in which the electron number density is 
not substantially altered by the ionization balance in the ele- 
ment E under consideration. 

It is useful to consider some special cases in relation 
to the structure of expression <32> . Consider the condition in 
which there is one dominant ion stage io, with Qi a « 1. The 
normalizations are such that as the stellar intensity is made 
either larger in amount or harder in energy, for a given base 
density no, the values of <5i will increase. As the dominant ion 
moves to higher stages, the lower ion stage fractions will be- 
come less significant. 

Consider the expression for the dominant stage, as given 

by 

S-i ■ ■ ■ Si -i 



Qi = 



(33) 



5- 1 +S- 1 S + ■■■ + (S-i ■ ■ ■ Si-i) 

Stages with i > io are not dominant, which can only mean that 
the photo-ionization rate is not large compared to the recom- 
bination rate, and so Si <C 1 for i > io - Since the other factors 
are large compared to unity (by definition, otherwise the stage 
io could not be dominant if lower stages were not trace), then 
in the denominator, the normalization must be dominated by 
the term that involves the product 5_i...5i _i, and so clearly 
Qi a ~ 1 since the numerator and denominator have like fac- 
tors. 

How then will adjacent stages scale with the density? For 
the next lower stage and using the above argument for the 
dominant term in the denominator, the ion fraction will be 

S-i 



oc n ( 



5~i ■ ■ ■ &o-i Si -i 
Similarly, the next higher stage will scale as 
S-i ■■■ S ia 



Qi -\ 



5-i ■■■ S io - 



5 io oc [n (9,)] 1 . 



(34) 



(35) 



Thus we arrive at the general rule, in the case where there is 
clearly a single dominant ion stage, that 



Qi oc [MO*)]* ' 



(36) 



This is not an especially fresh revelation. For example, the re- 
sult of equation I36i with its dependence on io — i for the 
case of spherical winds can be found in spirit in the discussion 
of Lamers, Cerruti-Sola, & Perinotto (1987). (The scaling ap- 
pears more clearly in Bjorkman et al. 1994 in their eq. [17].) 

A latitude-dependent ionization balance is important for 
the bi-stability effect (Pauldrach & Puis 1990) for generat- 
ing axisymmetric winds from the influence of stellar rotation 
(Lamers & Pauldrach 1991), although in this case it is a combi- 
nation of a latitude-dependent stellar radiation field, mass-loss 
rate, and terminal speed that generates the axisymmetric wind 
density, whereas we are concerned with a 'generic' latitude- 
dependent density and its feedback on the ionization balance 
in the wind. 

The radiative hydrodynamic models of Petrenz & Puis 
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(2000), who modeled axisymmetric line-driven winds from ro- 
tating stars, provide ion fractions with latitude (their Fig. 16). 
In their models the stellar radiation field is also latitude- 
dependent, and so impacts the ionization balance in addition to 
the varying wind density profile with latitude. Those authors 
are interested mainly in the inner wind region, where recombi- 
nation and resonance lines form, and they do not consider the 
formation of forbidden lines at large radius. 

Although we could have allowed for a stellar radiation 
field that was latitude-dependent, we have chosen to take the 
stellar radiation field as isotropic and focus on the large-scale 
asymptotic properties of the wind. But it should be noted 
that the latitude-dependence of the stellar radiation field and 
the latitude-dependent mass flux will not be independent, and 
could act to enhance or suppress the effects presented here, the 
question of which trend being dependent on the physics of the 
mass-loss processes and associated radiation transport. 
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Table 1. Atomic Data for Neon 



Recombination Rates Photo-Ionization Cross-sections d 
(lCT 12 cm 3 s" 1 ) (lCT 18 cm 2 ) 





<T = 58.9 


al = 0.283 x ^104/7; 


cr a = 14.6 


cQ = 1.71 x y/lO^/Te 


ct 2 = 2.7 


a% = 4.44 x y / lO i /T c 


as = 3.7 


a\ = 9.81 x ^/10 4 /T e 


0-4 = 0.54 


a\ = 25.6 x ^lOVTo 


cr 5 = 0.26 


<x\ = 42.7 x ^/10 4 /T C 


cr 6 = 0.49 


a\ = 65.8 x ^/l0 4 /T c 


err = 0.22 


ot% = 320 x ^/633/Tc 


cr 8 = 0.19 


a% = 616 x ^327/Tc 


o-g = 0.069 


ot c w = 1085 x ^183/T 





a From Osterbrock (1989) b Based on interpolation (see text) c From 
Verner & Ferland (1996) d From Verner et al. (1996) 



3.3 Model Results 

For computing a selection of model line profiles, we have cho- 
sen to focus on forbidden lines for the element of neon that are 
of particular interest for WR stars. The WR winds are selected 
because they show a number of strong forbidden emission 
lines in their IR spectra as revealed by ground-based and satel- 
lite observations (Barlow et al. 1988; Willis et al. 1997; Mor- 
ris et al. 2000; Dessart et al. 2000; Ignace et al. 2001; Morris, 
Crowther, & Houck 2004), and these are important for deter- 
mining metal abundances that can be used to test models of 
massive star evolution (e.g., Meynet & Maeder 2003). Neon 
is selected because lines from multiple stages of ionization of 
Ne can be observed in the IR, such as [Nell] and [Nelll], and 
possibly [NeV], 

Note that our models neglect clumping, which warrants 
a comment or two. Although clumping is well-known to af- 
fect diagnostics of WR winds, such as recombination lines 
(Hillier 1991) and radio emission (Abbott et al. 1981; Nugis, 
Crowther, & Willis 1998), and there is a growing body of evi- 
dence that clumping also influences the less dense O star winds 
(Crowther et al. 2002; Bouret, Lanz, & Hillier 2005; Fuller- 
ton, Massa, & Prinja 2006), a good theoretical understanding 
of wind clumping is lacking, because the numerical calcula- 
tions are rather complex (although there have been advances 
on this topic by Dessart & Owocki 2003, 2005). The details 
of how clumping initiates, where it initiates, how it evolves 
in the wind, and how it might vary in response to a latitu- 
dinal density distribution (not to mention the ionization bal- 
ance) remain more or less open questions. For spherical winds 
a constant clumping factor is straightforward to include, such 
as in Dessart et al. (2000) or Smith & Houck (2005). Assuming 
clumping did not vary with radius or latitude, it could be in- 
cluded in our expressions as a simple multiplicative correction 
factor, but we have chosen to ignore clumping in our models 
for now, because of the many questions currently surrounding 
the details of clumping. 



3.3.1 Ionization of Neon 

We have computed the ionization balance for Ne using the 
semi-analytic results of equation 1321 . TableQlists the atomic 
data used for these computations. The recombination coeffi- 
cients cii where the subscript i refers to recombination from 
Ne +t to the next lower ion stage, were taken from Osterbrock 
(1989) for i = 1 - 4. Values of qi for i = 8 - 10 come 
from Verner & Ferland (1996). The remaining coefficients 
were computed from an interpolation of these data. The photo- 
ionization cross-sections are values at the respective ionization 
edge frequencies and were taken from Verner et al. (1996). 
Here the subscript refers to photo-ionization from Ne +l to the 
next higher ion stage. 

We have simplified the calculations with two assump- 
tions: (a) the stellar radiation field is Planckian and (b) the 
photo-absorption cross-sections scales as o v = oi (vi/u) i , 
where o\ and z/j are values at the photo-absorption edge. Al- 
though these simplifying assumptions will change the detailed 
quantitative results of our models, they should not impact the 
qualitative results. For example, the scaling for the ion frac- 
tions presented in equation <36> should still hold. 

An example of the ionization balance is shown in Fig. [4] 
for neon plottted against the wind density scaling parameter 
n Cj o- The radiation field is taken to be Planckian with a temper- 
ature of 40,000 K, a reasonable effective temperature for WR 
stars derived from detailed spectral analyses (e.g., Hamann & 
Koesterke 1998). It is stressed that the ion fractions are plotted 
against n Cj o for different stars, and do not represent varying 
ion fractions as a function of wind density in a single star. For 
a given n Cj o, these are the ion fractions at every radius. The 
range of values for n e ,o are typical for WR winds, and the 
dominant ion varies from Ne IV at the lower densities to Ne 
III at the higher densities. At this temperature, Ne II is quite 
trace. Smith & Houck (2001) have found that [Ne II] 12.55 
microns appears in the spectra of some late WN and WC stars. 
For reference our models show that for a lower blackbody tem- 
perature of 30,000 K, Ne II becomes the dominant ion above 
2 x 10 13 cm 3 . 
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Table 2. Model Parameters 




11 11.5 12 12.5 13 13.5 

log n E0 (cm" 3 ) 



Figure 4. A plot of the ion fractions Qi for Ne+ l against the base 
wind density scale parameter n e o- Ne IV is dominant for the lower 
densities in this plot, but Ne III becomes dominant for n e o above 
about 3 X 10 11 cm -3 . As described in the text, these results assume 
a photoionizing spectrum from a 40,000 K blackbody. 

The results of Fig.|4]are admittedly a cheat. For the dense 
WR winds, changing the wind density scale impacts the emer- 
gent radiation field, which has feedback for the ionization 
balance in the wind. So it is not really correct to fix T e g at 
40,000 K and allow n e ,o to vary independently. We have com- 
puted ionization fractions based on WR spectral models pro- 
vided online 1 by the Potsdam group for WN and WC stars 
(Grafener, Koesterke, & Hamann 2002; Hamann & Grafener 
2003, 2004). These treat the radiation field and the wind den- 
sity profile consistently. For a broad sampling of the mod- 
els provided in terms of the transformed radius (Schmutz, 
Hamann, Wessolowski 1989) and stellar temperature, and us- 
ing equation I32i . the dominant ion stage is typically either 
Nell or Nelll, both being relatively strong forbidden lines 
commonly observed in the IR spectra of WR stars. Infre- 
quently, a higher ion stage is seen to be dominant, such as 
Nev. 

Although it would be more accurate for any given 
WR star to use the stellar radiation data from models like those 
of the Potsdam group (even better, the actual computed ioniza- 
tion data, but that information is not provided at their website), 
and to have better atomic data to relax assumption (b) above, 
ours is really an exploratory paper to investigate trends that 
axisymmetric winds will have for forbidden emission profiles. 
It is sufficient at this point to treat the latitude dependence of 
the density and ionization balance self-consistently within the 
stated assumptions to arrive at valid qualitative conclusions. 

For an axisymmetric wind, one has n c ,o a function of lat- 
itude. Fig. [4] indicates that changes in n c .o with latitude will 

1 www.astro.physik.uni-potsdam.de/PoWR.html 
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a Fig. 5 is for the equator denser than the pole; Fig. 6 is for the pole 
denser than the equator. See eqns. (24) and (25). 

lead to changes in the values of the ion fractions, possibly 
even in a shift of the dominant ion stage. So both the den- 
sity and the ion fractions will change. Generally, as the base 
density increases, the ion fraction for a stage i > io will tend 
to drop, thereby somewhat compensating for the effect of in- 
creased density on the emissivity. On the other hand if i < io, 
then increasing the base density will lead to a higher ion frac- 
tion, thus compounding the effect of an increased density for 
the emissivity. 

3.3.2 Line Profile Generation 

For an optically thin emission line formed over a large vol- 
ume, the total line luminosity from an axisymmetric wind will 
be independent of the viewing inclination, and given by the 
following integration: 

9 f°° 9 f +1 
L ; =8tt y rdr J j(r,6«) dfi«, (37) 

However, the emission profile shapes will depend on viewing 
inclination. 

We had previously defined stellar spherical coordinates. 
We now introduce stellar Cartesian coordinates (x*,y*, 2*) 
and observer coordinates (x,y,z). The observer is located 
along the positive z-axis. The viewing inclination angle i is 
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defined as the angle between z and the stellar symmetry axis 
z*, thus i — 0° is a pole-on view, and i = 90° is edge-on. We 
also choose y — y* without loss of generality. The observer 
has spherical angular coordinates (9, (f>). 

The emission profile shape is given by a volume- 
integrated emissivity for different isovelocity zones. The 
isovelocity zones are defined by 

Vz = -^oo COS 9 = —Voa fl, (38) 

and so are conical-shaped regions. Noting that dv 2 = 
—Vood/j,, combined with dv = i/q dv z /c to give du = 
(wtxj/Ao) dfi, an integral expression for the profile shape 
dLijdv becomes 

^> z )=4tt— [ r 2 dr f j(r,0,)#- (39) 
dv v °° Jr, Jo 

Evaluation of this integral must generally be carried out nu- 
merically (although see App.lAlfor the special pole-on case), 
using the spherical trigonometric relation that for 9 and <j> 
given for a point in the wind, 9* can be found from 

cos 9, = cos 9 cos i + sin 9 sin i cos (f>. (40) 

Clearly for an arbitrary viewing inclination i, a cone of fixed 
9 will cut across stellar latitudes from 9* = 9 — i to 9 + i, 
and so a particular isovelocity zone at 9 will sample different 
densities and ion fractions. 

3.3.3 Parameter Study 

A number of model forbidden line profile shapes have been 
computed for ions of neon, and these are shown in Fig.|5]for 
dense equatorial winds and Fig.|6|for dense polar winds. A key 
to the different models displayed there is provided in Table|2| 
The models assume a wind dominated by Hell at large radius 
with fi = 2, fj, e — 4, a total wind density scale of no 
fi x 10 12 cm -3 , and an effective blackbody temperature of 
T = 40, 000 K. Figure|4]shows that Ne III will be dominant, 
with Ne IV down by a factor of about 10, and other ion stages 
inferior in abundance by much larger factors. 

These parameters are roughly chosen to match the wind 
properties of WR 134 (Hamann & Koesterke 1998). The mo- 
tivation is that WR 134 is known to have an asymmetric wind 
(Underhill et al. 1990; Vreux et al. 1992; Schulte-Ladbeck 
et al. 1992; Harries, Hillier, & Howarth 1998). The emission 
line of [Caiv] 3.21 /j,m obtained from ISO and discussed in 
Ignace et al. (2001) hints at a double-peaked emission profile, 
although the signal-to-noise is not very good in this weak line. 
Unpublished spectra of the strong Hell 3.09 um recombina- 
tion emission line in WR 134 clearly shows a well-defined, 
albeit asymmetric, double-peaked shape, similar to emission 
lines shapes that have been reported by others. 

Although the basal wind density is consistent with the 
wind and star parameters reported by the Potsdam group, the 
stellar temperature of 40,000 K was chosen to approximate the 
correct ionization balance; in fact T„ reported by Hamann & 
Koesterke (1998) is about twice as high, but the Hell ioniza- 
tion edge causes the spectrum to deviate strongly from Planck- 
ian. As discussed previously, these are simply details in light 
of the primary goal of highlighting qualitative profile effects. 



There are six panels in Figs. [5] and [6] with each showing 
several line profiles. All of the emission profiles are plotted 
as normalized flux density against normalized Doppler shift. 
The flux density is related to the specific line luminosity via 
dLi/dv = 47rZ) 2 F v . For the normalization, the integrated flux 
density of the line is used so that each continuum-subtracted 
line-profile would have a flat-topped profile of unit height for 
a spherically symmetric wind. Note that each panel has a dif- 
ferent vertical scale, indicated at the right or left side. 

The panels illustrate the range of profile shapes that can 
result as different parameters are varied. Each panel is labeled 
(a)-(f) in the respective figures, and a key for the model pa- 
rameters with each profile is provided in Table|2| In that table 
the value of G for Fig.|5|is the ratio of the equatorial to polar 
densities at the same radius, but for Fig. |6| G is the polar to 
equatorial density contrast. In both figures, the top two panels 
are for the extreme viewing inclinations of pole-on (left) and 
equator-on (right), with a mild density enhancement of G — 3 
and a power exponent of m = 5. For the case of the dense 
polar wind seen at i = 0°, the profiles evolve from strongly 
double-peaked for Ne II (solid) to modestly centrally bubbled 
for Ne V (long dash). Equator-on, the trend is opposite. The 
parameters G and m are fixed, so these are examples of rever- 
sals in the qualitative shape of the profiles as feedback from the 
ionization dependence on the density scale with consequence 
for the line emissivity. For this pole-on case, an analytic so- 
lution to the profile shapes and these scaling effects can be 
derived, which is provided in AppendixlAl 

The two middle panels are for the fixed dominant ion state 
of Nelll, with m — 3 fixed, and G allowed to vary from 1 to 
9. The case G — 1 is spherical and so produces the flat-topped 
profile. The left middle panel is for a viewing inclination of 
i — 30°, whereas the right middle panel is for i = 60°. At 
lower left the exponent m is allowed to vary from 1 through 
8 as indicated, with G = 3 and i = 60° fixed for the Nelll 
ion. The panel at lower right overplots profiles of m = 3 and 
G — 3 for Nelll as the viewing inclination varies. 

There are several interesting comments and conclusions 
to be made from the parameter study: 

(i) The profile distortions from a flat-top shape for the 
dense polar wind models as compared to dense equatorial 
winds appear more severe for a given value of G and m; how- 
ever, one must interpret this with caution. Although at fixed m, 
both models have the same opening angle A9, this does not 
translate to the same solid angle. For example, if m — 1 for a 
polar wind model, then A9 = 60° (see eq. 1261 ). The relative 
solid angle of the dense polar wind component is AO/47T = 
0.5. For a model with a dense equatorial wind, a structure of 
the same solid angle using the density parametrization of equa- 
tion I24i would have A9 = 30° requiring m = 4.8. So for 
example to make a fair comparison of profiles between Figs. 
[5] and [6] as m is varied, it is the short dashed line of|5Je) that 
should be contrasted with the solid line of[6te). In so doing, 
one finds that indeed the profile shapes are different between 
models of dense equatorial and dense polar winds, but the dis- 
tortions from flat-top are comparable in amplitude. 

(ii) Panels (a) and (b) for both Figures shows that quite 
strong centrally peaked or double-homed profiles do result 
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Equator Denser Than Pole 




Figure S. Model forbidden line profiles for axisymmetric winds that are more dense at the equator than the pole. Table |2] provides a listing of 
the model parameters and ions of neon being modeled for each of the line types in the different panels. The ordinate is normalized flux density 
(such that all of the emission lines have the same equivalent widths), whereas the abscissa is normalized doppler shift. Note that the vertical scale 
is different for each panel as indicated. The axisymmetric wind is taken to have i> ma x the same for all latitudes, and u b s varies from —Vmax from 
the left side of a profile to +Vmnx at the right side of a profile. Each profile drops to zero line emission at the extreme velocity shift values. 



for the dominant ion (dotted line type in these panels) for 
favorable viewing perspectives of pole-on or equator-on. On 
the other hand, mid-latitude views as shown in panels (d) and 
(e) for i = 60°, statistically a probable viewing angle, reveal 
that distortions from flat-top shape are pretty mild, usually less 
than 20% deviations for the range of values in G and m shown. 
This is a result of the isovelocity zones, cones in this case, 
cutting across a fairly broad range of stellar latitudes and so 
sampling both high and low density regions. The end result is 
compensation such that the profile tends toward flat-top, unless 
the values of G or m become somewhat large. 

(iii) The profiles have been plotted in normalized form to 
emphasis the impact of model parameters on line shape. How- 
ever, both G and m also impact the total line flux. If the emis- 
sivity were just linear in density, then the total line emission 
would be conserved even though the profile shape were to de- 
viate from flat-top. The emissivity is only linear in density for 
the 'flat' part; it is quadratic in the 'linear' segment of the con- 



tribution function (see Fig.0. As a result, line flux is not con- 
served. Indeed, it may be greater or smaller (depending on the 
ion stage) than an equivalent spherical wind of the same total 
mass-loss. For the dominant ion, line flux will be increased be- 
cause of deviations from spherical flow. In general this biases 
ratios of line fluxes for deriving relative abundances. Specifi- 
cally, if a line ratio is formed from two ions of different species 
but with the same value of Ai (and Ai is constant with lati- 
tude for both ions), the effect cancels, and ion fraction ratios 
will be the same as if the wind were spherically symmetric. 
So valid abundance ratios can be derived even if the wind is 
axisymmetric, a consequence of the line being optically thin 
and stellar occultation effects negligible. However, this is only 
true if Ai is the same for both lines, otherwise the ion fraction 
ratio will deviate from what would have been obtained from 
an equivalent spherical wind. 

(iv) The model results clearly show to what extent line pro- 
files will deviate from flat-top as model parameters are var- 
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Figure 6. As in Fig.|6] but now for winds that are more dense at the poles than in the equator. Again, a key for the different line profiles is to be 
found in Tablel2l 



ied, but what can be concluded about the wind structure given 
an observed line profile? Without some a priori knowledge, 
there are qualitative degeneracies which combined with lim- 
ited signal-to-noise data can make the inverse problem chal- 
lenging. For example, changing ion stage (Ai) has effects that 
are similar to changing G and/or m, even toward making it 
ambiguous as to whether the wind is denser at the poles or the 
equator. Normally, based on the spectral properties, one does 
have an idea about the dominant ion stage. 

If the ion state is roughly known, then increasing G and m 
tends to make the profiles display stronger deviations from a 
flat-top shape. Both the total line flux and the amplitude of 
the features (central peak or double-horns) must be matched 
simultaneously as G and m vary. The width of these features 
must also be matched, which may discriminate between dense 
polar versus dense equatorial models. The viewing inclination 
also influences the sharpness of the features as well as their lo- 
cation within the profile. Additional data, such as nebular mor- 
phology, polarimetry, or possibly even rotational v sin i values, 
may help to limit the range of viewing inclination. In practice, 
if the viewing inclination is intermediate of pole-on and edge- 



on (as it likely will be), and the wind distortion is relatively 
mild, it may be hard or impossible to confidently determine 
the wind distortion parameters unless both the signal-to-noise 
and spectral resolution are quite high. 



4 DISCUSSION 

The objective of this study has been to consider how forbid- 
den emission line profiles might be used as a diagnostic of the 
large-scale structure of stellar winds that deviate from spheri- 
cally symmetric outflow. To this end, we have introduced para- 
metric expressions for dense polar or dense equatorial axisym- 
metric winds. The stellar radiation field and the wind terminal 
speed have both been treated as isotropic, but these are as- 
sumptions that should be relaxed in producing synthetic line 
profiles for a particular model wind. If the mass-loss from the 
star varies with latitude, it is to be expected that the terminal 
speed will also vary, possibly correlating or anti-correlating 
with the mass-loss rate depending on the specifics of the wind- 
driving physics. In combination with the effects of the ioniza- 
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tion balance, deviations of the line profile shape from flat-top 
may be suppressed or enhanced compared to the results pre- 
sented here. Relevant are the effects of clumping and velocity 
dispersions, such as those investigated by Runacres & Owocki 
(2002) for an O star model at large radius in the flow (100's of 
stellar radii). Exactly how such properties will vary will global 
asymmetries is unclear. However, for a spherical wind, if such 
effects persist over the scale r c , they should be detectable in 
the form of sloping, in contrast to more nearly vertical, high 
velocity wings to the forbidden line profile shape. For aspheri- 
cal winds, the expressions in this paper for thin forbidden lines 
can be modified to allow for a range of flow speeds and densi- 
ties in a "cell-by-cell" approach to doing the numerical volume 
integration for the emission profile. 
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A particular focus of this work has been to highlight the 
coupling between density and ion stage which are coupled in 
latitude. These can act through the line emissivity to exac- 
erbate or depress deviations of the line profile shape from a 
flat-top. Including the effects of a latitude-dependent radiation 
field and terminal speed may modify the conclusions that we 
have drawn from our assumed model, but the overall principle 
that resolved forbidden emission lines can be excellent tracers 
of deviations from spherical flow remains robust. At the same 
time, Figs.|5|and|6|do indicate that for winds with modest dis- 
tortions from spherical as viewed somewhat more equator-on 
(the case i = 60°) than pole-on (i = 30°), the line profiles 
of the dominant stage show modest deviations from a flat- 
top. Combined with instrumental smearing and finite signal- 
to-noise, such profiles might pass as essentially flat-top. 

As an observational strategy, it would be best to obtain 
forbidden emission line data for different ion stages within the 
same atom. This is not too common however, as ion stages sub- 
ordinate to the dominant one can have much smaller ion frac- 
tion values and produce only weak lines. The lines of [Nell] 
12.86 fim and [Nelll] 15.56 /im are an example of a pair of 
forbidden lines in adjacent ion stages that can both be seen in 
WR star winds (see App.lElfor a discussion of the applicability 
of the two-level atom approximation for Nelll). ISO observed 
WR stars with high spectral resolution at the appropriate wave- 
lengths, but only a few such stars were observed. One exam- 
ple showing both forbidden lines is WR 11 (van der Hucht 
et al. 1996), but this is a colliding wind binary system. Both 
lines were also seen in the luminous blue variable, P Cygni 
(Lamers et al. 1996). Already Spitzer has observed a number 
of WR stars with the IRS; unfortunately, the IRS generally 
has inadequate spectral resolution to analyze the line profile 
shapes to apply the results of this paper. An exception to this 
might be the extremely fast, but also extremely rare, winds 
of the oxygen-rich WO subtype ((Barlow & Hummer 1982; 
Kingsburgh et al. 1995). These winds sport terminal speeds 
upwards of 3000 km/s (e.g., Drew et al. 2004). Forbidden lines 
observed at high resolution with the IRS for the WO stars may 
be suitable to determine or place limits (subject to the signal- 
to-noise) on the geometry of the asymptotic structure of these 
winds. These are hotter stars, so instead of [Ne III], perhaps IR 
lines of [Ne V] will be present. 



APPENDIX A: ANALYTIC LINE PROFILES FOR A 
POLE-ON INCLINATION 

The latitude dependence of the axisymmetric density distribu- 
tion is described by 

/(ft.) = n(r, 8«)/n sph (r) = H' 1 [1 + (G - 1) g{8,)] , (Al) 

where H is a normalization set by eq. \21\ , and g(8 t ) = 
cos m 8, for a dense polar wind or sin m ft* for a dense equato- 
rial wind. 

In the case of a pole-on viewing inclination with i = 0°, 
one has that 8 = 6,. The emission line profile shape is then 
determined by 



dLi f . . 2 , 
oc / n2 (r)r ar 

dv Jr. 



n e Qi 2 , 
■ r ar 



i r, l+ n c/n c 
Using eq. 1361 and defining Ai = i — io, we have that 



dv 



iv , l + /r»/rg 
i 



dr 



f~ Ai (iS + f^uiy 1 du 



(A2) 
(A3) 

(A4) 
(A5) 



where the last line involves a substitution of u = R*/r and 
the lower limit to the integral assumes u c — R*/r c <C 1. This 
last integral has an analytic solution: 



dv 



lf(6* 



(A6) 



Since 8 = 8,, and the observed doppler shifted velocity in the 
profile is v z = v ma x cos#, one has that cos#» = ^z/v maJC , 
that can be substituted in for g{8,) to derive an analytic func- 
tion for the emission profile shape. 

Example profiles based on equation <A6> are shown in 
Fig. IA1I The upper panel is for a dense polar wind, and the 
lower one for a dense equatorial wind. Similar to Figs.|5]and 
|S| the line types are for Ai = — 1 (solid), (dotted), +1 (short 
dash), +2 (long dash). In those Figures, [Ne III] was the dom- 
inant ion stage for neon, but here Fig. IA1I is for any metal 
species. Unlike Figs. [5] and [6] a profile for Ai = —2 has been 
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Figure Al. Analytic solutions for forbidden line profile shapes as seen pole-on. Upper is for a dense polar wind, and lower is for a dense equatorial 
wind. For the dense polar wind of the upper panel, profiles progressively shift from a strongly double-horned emission shape with Ai = — 2 to one 
that is mildly centrally peaked at Ai = +2. The exact opposite progression occurs for a dense equatorial wind in the lower panel. Values of m = 5 
and G = 3 were adopted, and the profiles are normalized so that a spherical wind would produce a profile of unit height in this figure. 



computed for Fie. lAll and is shown as dot-dash line type. In 
all these cases, values of G — 3 and m = 5 were used. The 
profiles are normalized by their total emission. 

In the case of the dense polar wind as viewed pole-on, the 
profiles for the dominant ion stage and lower stages are sub- 
stantially double-peaked, because / is large along the poles. 
Even though G has a modest value, these profiles depend on 
/ in powers of 1.5, 2.5, and 3.5, so that strong double-horned 
profiles still result. In contrast, the profiles at Ai = +1 and 
+2 show a slight double-horned structure in the first case and 
a centrally bubbled one in the second. Both of these appear 
nearly flat-topped owing to the weak dependence on / at /~ ' 5 
and / +0,5 , respectively. 

In contrast, the profiles for the dense equatorial flows 
(lower panel) are all centrally peaked, except for one that 
shows a modest double-horned morphology. Here the denser 
flow is mostly in the plane of the sky corresponding to ob- 
served velocity shifts around line center. For the dominant ion 
stage and lower, this leads to modest or strong central peaks 



relative to a flat-topped profile and substantially depressed 
emission wings. Interestingly, the depressed wings could make 
the determination of the true line width problematic owing to 
uncertainties in the continuum placement for modest signal-to- 
noise data. Once again, as seen in the upper panel, the profiles 
for 1 and 2 stages above the dominant ion approach a some- 
what more flat-top appearance, and switches from centrally 
peaked emission for Ai = +1 to the double-horned appear- 
ance for Ai = +2. 



APPENDIX B: NEIII AND THE TWO-LEVEL ATOM 
APPROXIMATION 

The derivations of this paper are predicated on the assumption 
of a two-level atom. A forbidden line of considerable inter- 
est from WR stars is the [Nem] 15.55 ^m line ( 3 Pi - 3 P ), 
which is seen to be strong in a number of WR stars from ob- 
servations by ISO (Morris et al. 2000; Dessart et al. 2000) and 
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Figure Bl. A plot of normalized level population number densities 
for the five lowest energy levels in Nelll. Assuming highly NLTE con- 
ditions, the population in all levels greater than these five are taken to 
be negligible. The transition 2 — 1 corresponds to the forbidden line of 
[Nelll] 15.55 fim, and this figure shows in large part, this atom can be 
reasonably approximated as a two-level atom for conditions relevant 
to WR star winds. 



Spitzer (Morris et al. 2004; Cassinelli, private comm.). How- 
ever, Nelll has five low-lying energy levels, that are (in order 
of increasing energy): 3 Po, 3 Pi, 3 P2, 1 D2, and 1 So- The ques- 
tion is then, how bad is the two-level approximation for Nelll? 

Population number densities for these five levels, denoted 
1-5 in the order listed above, have been computed, and the re- 
sults are displayed in Fie. IB II To derive the populations, the 
equations of statistical equilibrium were solved allowing for 
collisional excitation and de-excitation, and radiative decay. 
Being interested primarily in the large radius conditions, the 
electron temperature was taken to be T c = 10, 000 K. Atomic 
data were taken from Pradhan & Peng (1994). Absorption and 
stimulated emission were included, assuming a blackbody ra- 
diation spectrum, but as expected, these were found to have no 
impact on the level populations. 

Fig. IB II shows that n\ + ni ~ n(Ne III) is reason- 
ably accurate, with an error not exceeding 10% at densities 
higher than the critical density, and for densities near and be- 
low the critical value, the error becomes negligibly small as 
n\ — > n(Ne III). This conclusion suggests that if Ne III is 
the dominant ion, then neon abundances can safely be derived 
from the 15.55 line in the two-level atom approximation, 
with a correction that is probably comparable or less than mea- 
surement errors associated with continuum placement for de- 
termining the line flux or errors in the distance measurement 
of the star or of other fundamental wind parameters (e.g., the 
mass-loss rate). 
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